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Capacitively coupled nano conductors 

Ratchet currents and exchange fluctuation relations 
Robert Hussein and Sigmund Kohler* 


We investigate electron transport in two quantum cir¬ 
cuits with mutual Coulomb interaction. The first circuit is 
a double quantum dot connected to two electron reser¬ 
voirs, while the second one is a quantum point contact 
in the weak tunneling limit. The coupling is such that an 
electron in the first circuit enhances the barrier of the 
point contact and, thus, reduces its conductivity. While 
such setups are frequently used as charge monitors, 
we focus on two different aspects. First, we derive trans¬ 
port coefficients which have recently been employed 
for testing generalized equilibrium conditions known as 
exchange fluctuation relations. These formally exact 
relations allows us to test the consistency of our mas¬ 
ter equation approach. Second, a biased point contact 
entails noise on the DQD and induces non-equilibrium 
phenomena such as a ratchet current. 


1 Introduction 

By now most quantum dots are designed with a nearby 
quantum point contact (QPC) whose conductance is af¬ 
fected by the charge state of the dot [1,2]. Then the QPC 
can act as monitor for the dot charge and provide time- 
resolved information from which the counting statistics of 
the electrons flowing through the quantum dot [3,4] and 
current correlation functions [5] can be reconstructed. 
When several quantum dots are strongly tunnel coupled, 
the wavefunction of their electrons becomes delocalized. 
Then charge detection corresponds to the measurement 
of the electron position and causes decoherence and lo¬ 
calization. A quantitative analysis [6] revealed that for a 
double quantum dot (DQD), good measurement correla¬ 
tions can be obtained only at the expense of a backaction 
strong enough to turn coherent inter-dot tunneling into 
classical hopping. 

Charge transport in a QPC in the weak tunnel limit con¬ 
sists of uncorrelated events in which an electron jumps 


from one lead to the other [7]. In technical terms, it rep¬ 
resents a Poisson process whose fluctuations are non- 
thermal and known as shot noise [8]. When acting upon a 
quantum system, these fluctuations not only cause deco¬ 
herence, but also excite electrons and drive the system out 
of equilibrium. In this way, the QPC may play a construc¬ 
tive role. In a asymmetric DQD (upper circuit in Fig. 1), 
shot noise may induce a ratchet current [9]. 

Coupled conductors can also be used for testing ex¬ 
change fluctuation relations [10-13] which are general- 



Figure 1 (a) Double quantum dot (DQD, circuit 1) coupled to 
a quantum point contact (QPC, circuit 2) operated in the weak 
tunneling regime. The QPC bias voltage V shifts the chemical 
potential of the right lead by p = -e 0 V. (b) Correlation function 
C(o0 = C + (ro) + C_M of the operator A by which the QPC 
couples to the left quantum dot for zero temperature (solid line) 
and k B T = 0.24eV (dashed), see Eqs. (4) and (8). 


* Corresponding author E-mail: sigmund.kohler@icmm.csic.es 
Instituto de Ciencia de Materiales de Madrid, CSIC, 28049 Madrid, 
Spain 


Copyright line will be provided by the publisher 


1 

































R. Hussein et al.: Capacitively coupled nano conductors 


ized equilibrium relations based on the assumption that 
each lead is in a Gibbs state. Then forward and backward 
rates are related by a Boltzmann factor which provides 
relations between transport coefficients. A prominent ex¬ 
ample is the Johnson-Nyquist relation between the con¬ 
ductance and the zero-frequency limit of the current cor¬ 
relation function [10,11]. While exchange fluctuation rela¬ 
tions are exact, practical computations of transport prop¬ 
erties often rely on approximations such as perturbation 
theory in the dot-lead tunneling, which may violate exact 
formal relations. In turn, exact fluctuation relations can 
be used to test the consistency of theoretical methods. In 
this spirit, it has been shown [14] that the Bloch-Redfield 
master equation [15]—despite being a successful, widely 
applied, and fairly reliable approach—is not always fully 
compatible with exchange fluctuation relations. 

In this work, we address two aspects of transport in 
capacitively coupled conductors: First, we compute a 
set of transport coefficients which can be used to exper¬ 
imentally test fluctuations relations along the lines of 
Refs. [16,17]. Moreover, we investigate to which extent 
our master equation results agree with exact formal rela¬ 
tions. Second, we characterize the non-equilibrium cur¬ 
rent through a DQD induced by the coupling to a QPC 
and study cross correlations between these subsystems. 


2 Model and master equation approach 

2.1 DQD coupled to a QPC in the tunnel limit 

The model sketched in Fig. 1 consists of two electric cir¬ 
cuits, the upper one being a DQD coupled to electron 
source and drain. We describe the DQD by two levels 
tunnel coupled with matrix element ft and detuned by e 
such that the one-electron states are split by A = Ve 2 + Q 2 . 
Thus, 

€ + + F2 + + 

Hdqd = ~{c' 2 c z - c[Ci) + —[C[C 2 + c' 2 Ci) + Un±n 2 , (1) 

where the fermionic operators c\ and c 2 annihilate an 
electron on the respective dot. The last term of F/dqd ac¬ 
counts for inter-dot Coulomb interaction with strength 
U, while we assume that the repulsion within one dot in¬ 
hibits double occupation. We will not study spin effects 
and, thus, work with spinless electrons. The leads are de¬ 
scribed by Hleads = Y. q eq{n Ll c, + n R) q ) with n a ,q = C f a ,qC a ,q 
the occupation number of mode q in lead a = L,R. The 
coupling between the leads and the dots is given by 

#DQD-leads = X C\ + £ VR,qc\ C 2 + h.C., (2) 

q ’ q 


where “h.c.” denotes the Hermitean conjugate. The ef¬ 
fective coupling is given by the spectral density T = 
2nY.q I VLi Ri q\ 2 8(e-e q ) which we assume energy indepen¬ 
dent. It provides the dot-lead rate TIh. 

Our second system is a QPC between two leads mod¬ 
elled by the Hamiltonian Hqpc = 'Lk e kc] c c k + 'Lk'£k'c' k ,c k ', 
where k and k 1 label the modes of the left and the right 
lead, respectively. The leads are weakly coupled by the 
tunnel Hamiltonian A = A + + A_, where 

A + — ^ Cfc, (3) 

k,k’ 

transfers an electron from the left to the right QPC lead, 
while A_ = A+ describes the opposite process. In a 
continuum limit, the matrix elements l k k : are encom¬ 
passed by the energy-independent QPC conductance 
G = 27rX J tfc/|t J t J fc/| 2 5(e-eQd(e-e)tO in units of the con¬ 
ductance quantum Go = e^/h. If an electron resides on 
dot 1, Coulomb repulsion enhances the barrier between 
the leads and, thus, reduces the tunnel matrix elements 
Lkk<- This effect is captured by a prefactor x = (1 - snj ) in 
the tunnel Hamiltonian such that 

Hqpc = *( A + + A-) = (1 - sn i) (A+ + A_) (4) 

accounts for both the QPC and its coupling to the DQD. 
For consistency, the dimensionless coupling s must obey 
s< 1. 


2.2 Master equation and full-counting statistics 

Our theoretical description is based on the formal elimi¬ 
nation of all four leads such that we remain with a reduced 
master for the DQD. After transforming the Liouville-von 
Neumann equation for the total density operator into the 
interaction picture with respect to Hdqd and the lead 
Hamiltonians, we derive a Markovian master equation 
[18] that captures the remaining terms to second order, 

OO 

p = -^[H s ,p]-— £ J dttii eads [V n ,[V n (-t), p ® R 0 ]], 

n o 

(5) 

where p is the reduced DQD density operator. Rq refers 
to the grand canonical ensemble of each lead a with a 
chemical potential shift p a measured with respect to a 
common Fermi energy e R . The operators V n represent 
H^ u p n c and the two tunnel contributions of H D QD-ieads- 
While the reduced master equation (5) fully describes 
the DQD, the lead degrees of freedom are traced out, so 
that information about the transported electrons gets lost. 
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To be able to recover it, we multiply before tracing out 
the leads the full density operator p ® Rq by a phase factor 
e' x ' N , where the elements of the vectors x = ( Xl>XR’ Xqpc) 
and N = ( Nr , Nr, IVqpc) refer to the left and the right lead 
of the DQD and to the QPC, respectively. 1 Then we ob¬ 
tain the master equation p = ££{%)p with the generalized 
Liouvillian S£{x)p = -(i//z)[H DQD ,p] +if DQD . leads (^)p + 
-^qpcOk) P- The now ^-dependent p obeys trp = < e' x ' N ), 
i.e., it is a generating function from which the moments 
of the lead electron distributions can be computed by 
taking derivatives with respect to components of X- Ac¬ 
cordingly, Z(x) = Jylntr p generates current cumulants, 
while derivatives at different times provide current corre¬ 
lation functions. After multiplication with the appropriate 
power of the electron charge -eo, we obtain the currents 
as the corresponding change of the lead electron number: 
h = ~ie 0 dZIdxQpc and /, = {-ie 0 l2){dZldxL~dZldxR) 
evaluated at x = 0. The definition of I\ is motivated by 
displacement currents in the leads which have the conse¬ 
quence that primarily the symmetrized current is experi¬ 
mentally accessible. While being irrelevant for the average 
current, this affects time-dependent quantities such as 
correlation functions [8]. 

The evaluation of i^DQD-ieads (%) for the incoherent 
DQD-lead tunneling is rather standard, see e.g. the ap¬ 
pendix of Ref. [19]. It essentially consists of jump oper¬ 
ators between many-particle DQD states that differ by 
one electron. The transition rates contain Fermi functions 
reflecting the initial occupation of the lead modes. The 
elimination of by contrast, is way less common 

and, moreover, describes the interaction of the two cir¬ 
cuits which is in the focus of the present article. Therefore 
it is worthwhile to discuss it in more detail. 

By evaluating the /-integral in Eq. (5) we end up with 
the QPC part master equation given by 

JZqpcW = ^qpc + (e^ QPC - D/ + + (e _ ' lQPC - 1)/“. (6) 

Its first term is the QPC Liouvillian 

1 r+°° 

-S^qpcP = 77 ^ / dtC{t)[x{-t)px + xpx{t) 

J— oo (7) 

-xx{- t)p - px{t)x\. 

The symmetrization of the time integral corresponds to 
neglecting energy renormalization stemming from princi¬ 


1 When an electron tunnels between the DQD and one of its 
leads, the DQD occupation and, hence, Nl + Nr change by +1. 
Therefore, a full description requires two independent counting 
variables. For the QPC, by contrast, /Vq PC + Mq PC = const, so 
that one counting variable is sufficient. 


pal values. A main ingredient to JZqpc is the correlation 
function of the QPC tunnel operator, C(t) = (A(/)A(0)) = 
C+ (I) + C_ (t), where C+ (I) = (A+ (f) A+ (0)) is readily evalu¬ 
ated from its definition and the assumption that the leads 
are voltage biased. In Fourier representation it reads [7] 


ho)±e oV 

(8) 

1 - exp [- (ha) ± eo V) / kg T] 

: T {-o))e [hM±eoV),kBT , 

(9) 


see sketch in Fig. 1(b). For energy absorption from the 
environment, it is evaluated at negative frequencies, for 
emission at positive frequencies. In the low-temperature 
limit, C(a») vanishes for hto < -\e 0 V\, which means that 
the size of the energy quanta absorbed by the DQD is 
limited by the QPC bias. Interestingly enough, C+ can 
be written in terms of e[l - exp (el kg T)] -1 , an expression 
appearing in master equations for Ohmic quantum dissi¬ 
pation. 

The second and third term in Eq. (6) stem from the 
counting field and vanish for Xqpc — 1 • 0. The superoperator 

1 f + °° 

J+p = -frzj dtC+{t)[x{-t)px + xpx{t)] (10) 

describes forward tunneling, while backward tunneling is 
given by the corresponding expression with C_ (/). 

For the numerical treatment, we decompose the su¬ 
peroperators into the eigenstates of Hdqd- In this basis, 
the interaction picture operators contain time-dependent 
phase factors, so that the time integrals yield delta func¬ 
tions and, finally, C± (to) has to be evaluated at the tran¬ 
sition frequencies of the DQD. Thus, we have obtained 
an explicit representation of the Liouvillian 5£ (%) which 
provides access to correlation functions [6], current cumu¬ 
lants [20], and transport coefficients [14]. Details for each 
method can be found in the respective quoted reference. 


3 Exchange fluctuation relations 

For typical system-lead models, one can write a formally 
exact expression for the cumulant generating function 
Z beyond the present perturbation theory. While gener¬ 
ally such expressions cannot be evaluated exactly, they 
provide formal properties such as the so-called exchange 
fluctuation relation [11] 

Z{x) = Z{-x-ililk B T ), (11) 

where Z(x) implicitly depends on the chemical poten¬ 
tials p. The derivation of this relation is based on the 
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assumption that each lead consists of a continuum of 
modes initially in a Gibbs state, while the central system 
(here: the DQD) possesses only a few degrees of freedom. 
A main interest in such relations stems from the fact that 
the Taylor coefficients of Z at p = x = 0 are transport co¬ 
efficients, i.e., experimentally accessible quantities such 
as the conductivity. Hence, Eq. (11) provides relations be¬ 
tween transport coefficients known as Onsager-Casimir 
relations [11]. 

To highlight the relation between the cumulant gener¬ 
ating function Z and Casimir-Onsager relations, we no¬ 
tice that by construction, the electron number in lead a 
changes by the current I a = -ie 0 dZ/dx a \x=o- Moreover, 
close to equilibrium p = 0, the current is linear in the volt¬ 
ages. Its slope is a transport coefficient, namely the (trans) 
conductance G a p = -d/ a /d/i^|/j=o which obviously is a 
second derivative of Z{%)- Taking the same derivative also 
on the r.h.s. of Eq. (11) yields S a pl kgT - G a p, where S a p 
is the (co) variance of the currents. Thus, we obtain the 
Johnson-Nyquist relation2 ksTG a p = S a p. Generalizing 
this concept, we define the transport coefficients 


K 


{-i) m e™ +n d m+n 
Pl " Pn ~~ dXaG" d Xa m dP(h ■ ■ ■ dpp n 


a v -a m 


Z(X) 


x=P=o 


( 12 ) 


by HdqdI<2> = E a \a). Thus, one often can work with the 
ansatz p a p = P a 8 a p which provides the Pauli-type mas¬ 
ter equation P a = Y. a ' ^aa'W^’a'- The transitions rates 
W aa i ix) follow from a basis decomposition of the master 
equation (5) after having introduced the counting vari¬ 
ables x = (Xl>XR’X qpc)- For our later reasoning, it is im¬ 
portant that the long-time solution of a Markovian mas¬ 
ter equation is dominated by the eigenvalue -X 0 (x) of 
the corresponding matrix W(x) that vanishes in the limit 
X — > 0 , i.e., the one that corresponds to the stationary 
solution in the absence of the counting variable. Then 
P ~ e~ /ln( T f which translates to Z = -Ao and implies that 
the generating function can be traced back to the com¬ 
putation of an eigenvalue of W[x ) [27]. Hence we can 
conclude that our RWA master equation complies with 
Eq. (11) if W(x) an d W(~x - ipl kg T) are isospectral. For 
the Liouvillian of the type ^DQD-ieads for the dot-lead tun¬ 
neling, this has already been demonstrated in Ref. [14]. 
Thus, we can restrict ourselves to the corresponding cal¬ 
culation for -SfQPcOf) in Eq. (6). 

By evaluating W aa Xx) f° r -S^qpc and,/+, see Eqs. (7) 
and (10), we find 

W aa Xx) = <? ! * QPC u^-a' + e~ lXQFC w~_ a " a? a! (13) 


Applying the same derivative to Eq. (11) provides general¬ 
ized Casimir-Onsager relations, which we abbreviate as 
Rcci- am _ q p Qr exam pi e = o denotes the Johnson- 

Pl—Pn r p 

Nyquist relation. 

Since approximation schemes are not necessarily con¬ 
sistent with exact formal properties, the question arises 
whether our master equation approach complies with 
Eq. (11). Its consistency has been verified for various spe¬ 
cific situations [21-26], while a general proof has been 
given only in the classical limit [12,13] and within rotating- 
wave approximation (RWA) in the single-particle limit [12] 
and for many-particle states [14]. Beyond RWA, Eq. (11) 
may be compromised. The deviations, however, tend to 
be tiny, in particular in the low-temperature limit [14] 
which is non-trivial since the temperature appears in the 
denominator. Here we explore the situation for a master 
equation with the non-conventional Liouvillian JZqpc- In 
doing so, we extend the results of Refs. [17,23] to the pres¬ 
ence of quantum coherence and those of Ref. [14] to the 
coupling to a QPC. 

3.1 Bloch-Redfield equation in RWA 

fn the limit of very weak coupling between the central 
system and the leads, the reduced density operator be¬ 
comes eventually diagonal in the energy basis defined 


with the QPC forward and backward rates 

K~a' = c ±( E a'-E a )\(a\x\<x')\ 2 . (14) 

Since the Liouvillian 5£{0) conserves the trace of the den¬ 
sity matrix, we can obtain the diagonal elements W aa 
from the normalization condition J_ a P a = 0 yielding 
w aa = -Za'^ia Wa'a(0). Using Eq. (9) for p = -e 0 V, we 
find for the QPC contribution by straightforward algebra 
the relation 

Waa' (-*qpc -iplk B D = e E ^' lkBT W a i a {xQPc)e~ EalkBT ■ 

(15) 

Together with the corresponding relation for the dot-lead 
tunneling [14], we obtain 

W(x) T = S~ l W{-x~iplk B T)S, s=e~ HD Q DlkBT , (16) 

where T denotes matrix transposition. This means that 
the RWA Liouvillian with the modified counting variable 
relates to the original one by a transposition and a sim¬ 
ilarity transformation which both leave the eigenvalues 
invariant. Therefore, we can conclude that the RWA mas¬ 
ter equation agrees with the exchange fluctuation rela¬ 
tion (11). 
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3.2 Deviations and experimental tests 

Despite possessing the desirable consistency with Eq. (11), 
a RWA master equation has obvious limitations when off- 
diagonal density matrix elements play a role. This is, e.g., 
the case for a strongly biased, but undetuned DQD in 
which the inter-dot tunneling is weaker than the dot-lead 
tunneling. Then an electron that enters from the source 
will stay on dot 1 for a while until it proceeds to dot 2. From 
there it will tunnel rapidly to the drain. Thus, the DQD will 
will predominantly be in the state cj|0), so that its den- o 
sity operator in the energy basis has a large off-diagonal 11 
contribution and the RWA will fail [28]. As a drawback, ^ 
however, a treatment beyond RWA is generally not fully § 
consistent with exchange fluctuation relations [14], which ^ 
motivates a quantitative study of possible discrepancies. -2 

Experimentally, the cumulant generating function is 
not directly accessible and, thus, one may test instead •% 
Casimir-Onsager relations derived from it. We here focus 
on transport coefficients that have been tested with an 
Aharonov-Bohm ring [16] and a DQD coupled to a QPC 
[17]. We start with the Johnson-Nyquist relation and find 
in accordance with Ref. [14] that it is fulfilled within our 
numerical precision. 

Thus, for revealing possible discrepancies, we have to 
employ Casimir-Onsager relations of higher order such 
as kg 'l'K ] u + K^ 1 = 0 which follows directly from R l u = 0 
and A’ 111 = 0. Figure 2(a) depicts the master equation re¬ 
sult for both terms involved in this relation. Their main 
contribution appears at the conductance peaks, i.e., at 
Fermi energies that obey H {n+V> = + ep, where H ln> is 

the energy of the lowest state with n electrons on the DQD. 
Here, E (01 = 0, £ tl) = -\Ve 2 + D. 2 , and £ (2) = U, as can be 
readily verified by diagonalizing Hdqd in Fock basis. The 
interpretation of this condition is that an electron from 
the Fermi surface can tunnel between the lead and the 
DQD without energy cost, so that already a small voltage 
can induce a current. With the scale chosen, kg T and 
/Q 11 appear identical as they ideally should be. A closer 
look, however, reveals a relative difference of the order 
1CT 3 . Such tiny differences are beyond the resolution of 
related experiments [16,17], which lets us conclude that 
our master equation treatment is sufficiently precise to 
compute the quantities employed for typical tests of ex¬ 
change fluctuation relations. 

Beyond such practical issues, it is interesting to see 
how these discrepancies scale with the system parame¬ 
ters. For the “conventional” dot-lead jump operators such 
as those in ii?DQD-ieads> one generically finds for small 
dot-lead rate R oc T 3 , while for small inter-dot tunnel¬ 
ing, Roc Q 2 [14]. For SA)pc, see Eq. (6), the relevant pa¬ 
rameters are the dimensionless coupling s and the QPC 




G 



(circles) computed with our Bloch-Redfield master equation 
as a function of the Fermi energy for coupling 5=1, QPC 
conductance G = 1, temperature T = OACl/kg, and the DQD 
parameters T = e = 2U = Q. While according to the exchange 
fluctuation relation must be differ only by a factor k B T , the 
master equation predicts an additional tiny spurious deviation 
(dashed line). The vertical lines mark the location of the con¬ 
ductance peaks as discussed in the text, (b) Scaling of the 
deviations for kgTK ^ + kfj 11 and kT 111 from their exact theo¬ 
retical value zero as a function of the QPC conductance for 
e P = U+M2. 


conductance G, which brings the scaling as a function 
of 5 and G to our attention. Figure 2(b) depicts the peak 
values of the deviation shown in Fig. 2(a). It indicates that 
kgTK + kfj 11 oc 5 2 G. In particular for s = 0.1, the l.h.s. of 
this relation involves computing the tiny difference be¬ 
tween two much larger numbers, which is at the limit 
of our numerical precision. Indeed one notices that for 
small values of G, the data is compromised by numeri¬ 
cal errors. We therefore also studied the Casimir-Onsager 
relation K ] 11 = 0 which follows from A 111 = 0. Since this 
relation consists of only one transport coefficient, it does 
not suffer from the mentioned numerical problem. The 
corresponding data in Fig. 2(b) confirms for the deviations 
the scaling oc G and oc 5 2 . 
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Thus we can conclude that the full Bloch-Redfield mas¬ 
ter equation for the QPC coupling is consistent with ex¬ 
change fluctuation relations up to corrections Roc s 2 G, 
i.e, corrections linear in G. By contrast, the mentioned 
scaling with the DQD parameters T and Q is more favor¬ 
able [14]. Nevertheless, the Bloch-Redfield approach can 
be safely applied for the computation of transport coeffi¬ 
cients up to second order and when a RWA treatment is 
sufficient. 


4 Tunnel contact as driving source 

When non-equilibrium noise acts upon the DQD, it will 
induce electron transitions from the ground state of the 
DQD to the excited state, see Fig. 1(a). Once in the excited 
state, the electron will leave predominantly to the right 
lead, while subsequently an electron from the left lead 
enters. In this way, the QPC can spurs a directed current 
in an unbiased DQD circuit. Its direction depends on the 
sign of e, which implies a current reversal point at e = 
0. For strong DQD detuning, this mechanism has been 
proposed as noise detector [29]. Going beyond that work, 
our master equation description starts from a complete 
model that includes the noise source. Moreover, it allows 
us to study also the backaction on the QPC as well as cross 
correlations between the circuits. 


4.1 Ratchet current 


To estimate the ratchet current, we compute the rates 
for the scenario sketched above. A main role is played 
by electron transitions from the one-electron ground 
state of the DQD, |g) = 11> cos0 + |2) sirup to the ex¬ 
cited state |e) = -|1) sin</>+ |2)cos (f>, where cos(2(/>) = el A 
and | n) = 10). These transitions are induced by the 

ffamiltonian (4), and occur with the golden-rule rate 
y = (s 2 /4) sin 2 (2(/))C(-A/h)/h. The first factor of this ex¬ 
pression stems from the matrix element (e| rq |g) and ac¬ 
counts for the delocalization of the DQD orbitals. 

When the electron has reached the excited state, it 
can leave the DQD through either lead, with a probabil¬ 
ity that depends on the overlap of | e) with the localized 
states |1) and |2). The corresponding probabilities read 
|(l|e)| 2 = sin 2 ([> and |<2|e>| 2 cos 2 0. Since these transitions 
contribute to the current with opposite sign, we finally 
obtain I = eoy(cos 2 <fi - sin 2 <p), which can be written as 


/ = 


eo s 2 e0 2 
h 4(e 2 + Q 2 ) 3/2 


CC-A/ft). 


(17) 



e [Q] 


Figure 3 Ratchet current through the DQD as a function of 
the detuning e for tunneling matrix element O = 10T, QPC con¬ 
ductance G = 1 and voltage V = 2.5Q/e 0 , coupling strength 
s = 0.1, interaction U = 500, and various temperatures. The 
dotted line marks the analytical solution (17) for kgT = 0. The 
vertical lines marks the condition e 2 + 0 2 = (eoV) 2 beyond 
which the DQD splitting is always larger than eoV. Inset: current 
maximum as a function of the QPC bias voltage V for kg T = 0 
(circles) and kgT = 1.20 (triangles) in units of eoQIH. 


This expression generalizes the result of Ref. [29] to the 
regime |e| < O, in which we expect current reversals as a 
most relevant feature for possible applications. 

Our numerical solution allows refining the picture 
drawn by the golden-rule calculation. Figure 3 shows 
the ratchet current as a function of the DQD detuning, 
which for zero temperature by and large confirms the 
behavior predicted by Eq. (17), but indicates that the an¬ 
alytic approach overestimates the ratchet effect. When 
the QPC temperature is increased, its current acquires a 
thermal component so that the effective driving becomes 
stronger. Accordingly, we witness an enhancement of the 
the ratchet current. By contrast, when both the QPC tem¬ 
perature and the DQD temperature are increased simul¬ 
taneously, the bias voltage plays a smaller role and both 
circuits are dominated by thermal noise. Eventually the 
system reaches thermal equilibrium in which the current 
vanishes. The dashed line in Fig. 3 shows that already for 
the moderate temperature T = 1.20/ ks, the ratchet cur¬ 
rent is significantly smaller than at T = 0. 

The limitations of our analytic approach are also vis¬ 
ible in the current maxima. While Eq. (17) for e 0 V » e 
predicts a linear increase with V, the data shown in the 
inset of Fig. 3 demonstrates a sub-linear growth. This de¬ 
viation is even more clearly visible in the ratchet current 
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oj [n/h\ 

Figure 4 Absolute value of the correlation coefficient r \2 be¬ 
tween the DQD current and the QPC current as a function of 
the frequency for the detuning e = Q. All other parameters are 
as in Fig. 3. The vertical lines mark the values o) = T/h and 
a) = Alh, respectively. The inset shows the ratchet current in 
dependence on the applied bias voltage for zero temperature 
(solid line) and for k B T = 1.2Q (dotted). 


as a function of the QPC bias V for constant detuning 
(inset of Fig. 4). The main reason for this discrepancy is 
that our analytic approach is based on the assumption 
that the intra-dot excitation is the slowest process and, 
thus, governs the dynamics. Upon increasing V, however, 
this assumption will be violated at some stage and we 
would have to employ rate equation description along the 
lines of Ref. [30] which considers also the de-excitation 
\e) — |g>. 

4.2 Current noise and cross correlations 

To characterize current fluctuations, we consider the 
correlation functions S a p{t- t') = ( I a (t)Ip{.t ') - (I a )(Ip), 
where a and /3 label the two circuits. For the auto¬ 
correlations we find in the frequency domain roughly 
Saah'd ~ (la) (not shown), i.e., the typical frequency- 
independent value for the shot noise of a Poisson process 
[31]. The physical picture behind this behavior is that of 
uncorrelated tunnel events which holds for a tunnel con¬ 
tact in the low-temperature regime as well as for the DQD 
in the regime in which the transitions |g) —* |e) represent 
the bottleneck for the electron flow. 

Let us turn to the cross correlation S 12 which contains 
information about the mutual influence of the circuits 
on each other. A qualitative picture of the physical pro¬ 


cess emerges from normalized correlations. In the spirit 
of full-counting statistics of mesoscopic transport, it is 
common to define the Fano factors F aa {/U)) = S aa (oi)// a 
and the cross Fano factor Futa)) = SuicujlVhh- These 
quantities relate to the g (2) function in optics and hint on 
bunching and anti-bunching of transport events [8,32]. 
In a typical experimental realization, the currents in the 
two circuits may differ by several orders of magnitude 
[9]. Therefore we are not so much interested in corre¬ 
lated counting, but rather in the correlation coefficient 
U 2 (<^) = S 12 (m) / \/Si 1 (w) S 22 (<t>) which is normalized to 
the frequency-dependent auto-correlations. In our case, 
ri 2 (at) and F 12 (oi) are practically the same owing to the 
prevailing Poissonian nature of each current and the cor¬ 
responding almost constant S aa (w). 

Figure 4 shows the correlation coefficient for various 
temperatures. It is characterized by two energies. Below 
the dot-lead coupling [ , we find a plateau after which a de¬ 
cay oc w” 1 sets in. At the level splitting of the DQD, A, we 
find a peak with a value up to r 12 (A/fi) » 0.85 indicating 
strong correlations between the ratchet current and the 
QPC at resonance. Thus, correlations mainly appear when 
the QPC can induce resonant transitions between DQD 
orbitals. The sharp Lorentzian form of the peak indicates 
the relevance of quantum coherence, which implies that 
this feature is beyond the treatment with rate-equations 
that leads to the analytic solution (17). For smaller values 
of A, coherence eventually gets lost and the peak sub¬ 
merges below the low-frequency plateau. As expected for 
a coherence effect, the peak becomes smaller at higher 
temperatures. 


5 Conclusions 

We have investigated a setup composed of two capaci- 
tively interacting nano conductors, namely a QPC in the 
weak tunneling limit and a DQD. While in most recent 
experimental realization of such setups, the QPC serves 
as charge monitor for the DQD, we were interested in 
aspects beyond that. In particular, we focused on the non¬ 
equilibrium DQD dynamics induced by the shot noise of 
the QPC current and on testing the validity of exchange 
fluctuation relations. 

For the theoretical treatment, we have derived a Bloch- 
Redfield master equation augmented by a counting vari¬ 
able for the electron number in each lead. This allows 
one to eliminate the leads within second-order perturba¬ 
tion theory without losing information about the trans¬ 
ported electrons. For the dot-lead tunneling, one obtains 
“conventional” jump operators, while the coupling to the 
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QPC yields dissipative terms resembling those stemming 
from the coupling to a bosonic heat bath. We demon¬ 
strated that within RWA, the resulting unconventional 
jump operators are fully consistent with exchange fluctu¬ 
ation relations. Beyond RWA, we numerically found de¬ 
viations which decay linearly with the QPC conductance 
and quadratically with the coupling. In the experimentally 
relevant regime, these deviations are rather small and do 
not inhibit computing with a Bloch-Redfield formalism 
transport coefficients for testing Casimir-Onsager rela¬ 
tions. The Johnson-Nyquist relation even turned out to be 
fulfilled within our numerical precision. 

When using the QPC as a driving source, one can in¬ 
duce a ratchet current in the DQD. As a feature most rel¬ 
evant for applications, it possesses a current reversal as 
a function of the detuning. While generally the ratchet 
current reflects the spectrum of the effective noise that 
the QPC entails on the DQD, our numerical results re¬ 
vealed significant deviations from the behavior found 
with a simple golden-rule treatment. The correlations be¬ 
tween the currents in both circuits are most pronounced 
at frequencies that correspond to the energy splitting of 
the DQD, which emphasizes the role of quantum coher¬ 
ence. For a typical inter-dot tunnel coupling of O = 25/ieV, 
the ratchet current is of the order of several pA, a value 
that can be measured straightforwardly with present tech¬ 
niques. Given that the setup investigated is common in 
today’s quantum dot design, we are convinced that our 
results will inspire future experiments. 
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